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ABSTRACT 

We present an analysis of CO emission lines from a sample of T Tauri, Herbig Ae/Be, and transitional 
disks with known inclinations, in order to study the structure of inner disk molecular gas. We calculate 
CO inner radii by fitting line profiles with a simple parameterized model. We find that, for optically 
thick disks, CO inner radii are strongly correlated with the total system luminosity (stellar plus 
accretion), and consistent with the dust sublimation radius. Transitional disk inner radii show the 
same trend with luminosity, but are systematically larger. Using rotation diagram fits, we derive, for 
classical T Tauri disks, emitting areas consistent with a ring of width ^0.15 AU located at the CO 
inner radius; emitting areas for transitional disks are systematically smaller. We also measure lower 
rotational temperatures for transitional disks, and disks around Herbig Ae/Be stars, than for those 
around T Tauri stars. Finally, we find that rotational temperatures are similar to, or slightly lower 
than, the expected temperature of blackbody grains located at the CO inner radius, in contrast to 
expectations of thermal decoupling between gas and dust. 
Subject headings: protoplanetary disks; stars: pre-main-sequence 



1. INTRODUCTION 

Observations of the inner regions (< 5 AU) of circum- 
stellar disks are essential for obtaining a complete un- 
derstanding of star and planet formation. The inner disk 
interacts with the central star, thereby controlling the ac- 
cretion and ejection of material and setting the timescale 
for star formation and disk evolution. In addition, inner 
disks are the likely birthplace of terrestrial planets. How- 
ever, these regions are difficult to study, because of their 
small angular size at the distance of nearby star-forming 
regions and the proximity to their parent stars. 

Recently, significant advancements in our understand- 
ing of inner disks have been made with two complemen- 
tary techniques: IR interferometry and high- resolution 
spectroscopy. IR interferometry is capable of observ- 
ing thermal emission from the dusty component of disks 
with resolutions of a few mas (probing size scales down 
to a few hundredths of an AU; Millan-Gabet et al. 
1999; Eisner et al. 2003; Akeson et al. 2005b, Mon- 
nier et al. 2005). These observations have shown that 
standard, opti cally thick accret ion disk models (such 
as outlined by iHillenbrandl I1992D fail to simultaneously 
fit spectral energy distributions (SEDs) and v isibilities 
(jMillan-Gabet et al.l [200TI : lAkeson et all I200I . Better 
fits to visibilities were obtained with an optically thin 
inner gas disk interior to a hot, ptiffed-u p wall of dust 
pullemond et al.ll2001l:lNatta et al.ll2001[ ). with the dust 
wall located where temperatures are high enough to sub- 



limate silicate grains onnier et al.ll2005D . Recent ob- 
servations and modeling efforts have introduced com- 
plications to this simple picture, including a signifi- 
cant source of emission of unknown origin inside the 
dust sublimation radius (fo r an extensive review, see 
iDullemond fc Monnieii[2010f) . However, the major con- 
clusion — that there is a defined inner radius consistent 
with dust sublimation — still holds true. 

The complementary technique of high-resolution spec- 
troscopy has been used to study rovibrational emis- 
sion from hot gaseous molecules , specifically H2O 
and CO ( e.g.. iBrittain et al.l 120031: iNaiita et all 120031: 
Blake & BoogertI 120041: iCarr et all 120041: iRettig et al.l 
120041: ISalvk et al.l l2008^. that originate in the inner disk 
atmosphere. Because of the high temperatures (> 1000 
K) required to populate vibrationally excited states, CO 
vibrational emission originates at disk radii similar to 
those probed by IR interferometry (< a few AU) . When 
the spectral resolution is high enough that the emission 
lines are spectrally resolved, the line profile acts as a 
proxy for the spatial location of the gas, under the as- 
sumption that the emission originates in a rotating Ke- 
plerian disk. Since the emission line wings represent the 
highest observed velocities, the flux in the line wings orig- 
inates at the inner edge of the molecular disk. Thus, 
emission lines (after a correction for disk inclination) can 
be used to measure molecular gas inner radii. 

In this work, we present observations of CO v—l^Q 
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emission from a large sample of protoplanetary disks, in- 
cluding "classical" optically thick disks around T Tauri 
stars (cTTs's) and Herbig Ae/Be (HAeBe) stars, as well 
as transitional disks — disks who se inner regions are 
depleted o f smaU dust grains fe.g.. ' Koerner et al1ll993l : 
iCalvet eT"al . 2002). The disks in the chosen sample all 
have known inclinations, allowing us to use line profiles 
to derive CO inner radii, which we then compare to the 
interferometric dust inner radii. We also fit the emission 
line fluxes with a local thermodynamic equilibrum (LTE) 
slab emission model, and relate the bulk properties of the 
gas to the measured CO inner radii. 

2. OBSERVATIONS AND REDUCTION 

Most of our spect ra were obtained with NIRSPEC 
(jMcLean et al.l Il998f ). a high-resolution {R -25,000, 
FWHM~12.5 km s"'^) spectrometer on the Keck II 
telescope. These observations are derived from sev- 
eral observing runs spanning the years 200I-20I0, and 
are part of a large NIRSPEC survey of protoplanetary 
disks, portions of which have been previously presente d 
(|Blake fc Booeertl [200l ISalvk et al.l [20071 [200l [20Tll) . 
The data were observed in the M band in echelle mode 
with a 0."43x24" slit. Each source was observed in at 
least two grating settings, thereby encompassing wave- 
lengths between — 4.65 and 5.15 /xm, with the exception 
of a hole between orders at ~ 4.8-4.95 fj,m. This wave- 
length range covers the first two R-branch lines and the 
low/mid P-branch (J = 1 - 12 and J = 30 - 40) of 
the w=l— i>0 CO rovibrational spectrum, as well as the H 
I Pf/3 and Hue transitions, and in a few cases the HnS 
transition. An observation log is provided in Table [T] 

TW Hya's CO emission lines were not well resolved 
with NIRSPEC and this source was therefore observed 
with Phoenix (Hinklc ct al. 2 003 ) on Gemini South, as 
first reported in ISalvk et al.l (|2007l ). It was observed on 
2006 April 7 and 8 using the 0."35x 14" slit. With its sig- 
nificantly higher spectral resolution {R ^ 60,000, FWHM 
~ 5 km s~^), but smaller spectral coverage, the Phoenix 
observations resolve the emission lines from P(6) to P(9). 

Objects were observed in nod pairs and then differ- 
enced. Exposure times (integration time multiplied by 
co-adds) were limited to one minute in length to min- 
imize atmospheric changes between nods. Wavelengths 
were calibrated using telluric emission lines. Nearby A 
and B standard stars were observed to correct for tel- 
luric absorption features. Any H I lines present in the 
standard star spectra were fitted by Kurucz models be- 
fore dividing source spectra by the standard. Standard 
stars were also utilized for flux calibration, using M- 
band fluxes estimated from Two Micron A ll Sky Survey 
X-band photometry (Skrut skie et al.l[2006l ) and spectral 
types from the literature. Portions of the spectra with 
poor atmospheric transmission (typically < 70%, but the 
exact percentage for each source was determined empiri- 
cally) were removed. Wavelengths were shifted to correct 
for a Doppler shift due to Earth's motion, which depends 
on the time of year. Many sources were observed at mul- 
tiple Doppler shifts, thereby "filling in" regions of poor 
atmospheric transmission and creating complete line pro- 
files, even for low-excitation lines. (High excitation CO 
lines do not suffer greatly from telluric absorption). 

A more detailed explanation of the data acquisition 
and reduction used for the NIRSPEC observations in our 



survey can be found in ISalvk et all (l2009l) . 

3. SOURCES AND SPECTRA 

The sample analyzed in this work consists of 32 young 
stars with circumstellar disks, including cTTs's, HAeBe 
stars, and stars with transitional disks (31 from our own 
survey, and one obtained from Najita et al. 2003). A 
description of the sample can be found in Table [2] This 
particular subset of our large survey was chosen for anal- 
ysis based on the following characteristics — CO in emis- 
sion, with no large line-shape asymmetries, and known 
disk inclination — with the goal of measuring CO inner 
radii from the line profiles. (Note that only two stars 
with disks — SR 24 and VSSG 1 — were eliminated due 
to line asymmetries, so such asymmetries are not charac- 
teristic of typical disks) . Our sample spans a large range 
of stellar masses and more than two orders of magnitude 
in luminosity, making this the first study capable of ex- 
ploring and comparing inner gas radii for a large range 
of luminosities and stellar types. In addition, our sample 
includes nine transitional disks (DoAr 44, CM Aur, HD 
135344 B, HD 141569 A, LkCa 15, LkHa 330, SR 21, TW 
Hya and UX Tau A) , making it possible to systematically 
compare their inner radii with classical disks. 

The complete set of spectra is shown in Figure [U 
with CO u = 1 — emission lines marked with dot- 
ted lines. Note that every source shows CO P- and R- 
branch emission lines (and most also show between one 
and three H I emission lines), but the spectra display a 
variety of line / continuum ratios and excitation tempera- 
tures (as reflected in the ratio of low- to high-excitation 
line strengths). 

4. ANALYSIS OF LINE PROFILES 
4.1. Constructing Line Composites 

To increase signal-to-noise ratio (S/N) and obtain good 
velocity coverage for line profiles, we constructed and 
analyzed line-profile composites. Because some of the 
lower-excitation CO lines can be contaminated by fore- 
ground absorption (e.g., see HL Tau in Figure[T|), and be- 
cause lower-excitation lines are more affected by telluric 
absorption, we created the line composites from higher- 
excitation ( Jup >20) lines only, whenever possible. The 
use of high-excitation lines also has the advantage that, 
since they originate from only the hottest gas, the inner 
edge of the disk makes a relatively large contribution to 
the total line fiux. However, when high-excitation lines 
were not available, we used others. There is sometimes, 
but not always, an increase in line width with excitation, 
which could potentially introduce biases into the analy- 
sis; we will discuss this in detail in Section 14.2.31 The 
set of lines used to construct each composite is listed in 
Table [3l along with the measured line width. 

The procedure for creating line composites was as fol- 
lows. Observed lines were screened for contamination 
(from ^■^CO and u = 2 — !> emission), and contaminated 
lines were eliminated. Remaining lines were centered at 
the theoretical line center (Rothnian ct al. 1992), inter- 
polated on a 5 (3 for TW Hya) km s^^ grid, and then 
averaged. Line composites are plotted in Figure O 

About half of the line composites are single-peaked, 
while half have evidence for somewhat double-peaked 
profiles. Two source composites, AA Tau and SU Aur, 
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Fig. 1. — Complete set of NIRSPEC spectra. Dashed vertical lines mark CO P- and R-branch lines. 4.654, 4.673, and 5.129 /im emission 
lines are HI Pf/3, Hue, and Hu5, respectively. Additional features present in some spectra include ^^CO 1^0 and I^CO 2^1 lines. 
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TABLE 1 
Observation Log 



Name 


Dates 


AA Tau 


2003 


Nov 3, 2004 Dec 27, 2004 Dec 30 


AB Aur 


2001 


Jan 30, 2001 Aug 7-8, 2002 Jan 3, 2002 Dec 18 


AS 205 N 


2002 


Apr 21, 2002 Jul 22 


BP Tau 


2009 


Oct 11 


DF Tau 


2005 


Dec 21 


DL Tau 


2003 


Nov 3, 2007 Dec 26 


DoAr 44 


2002 


Jul 22, 2004 Jul 23-24, 2005 Apr 26, 2010 Apr 22 


DO Tau 


2005 


Dec 18 


DR Tau 


2002 


Dec 17-18, 2005 Dec 18 


GG Tau 


2002 


Jan 3, 2002 Dec 17-18, 2003 Nov 2-3 


GK Tau 


2009 


Dec 27-28 


GM Aur 


2004 


Dec 27 


GSS 39 


2009 


Jul 13, 2010 Apr 22 


HD 135344 B 


2005 


Apr 24, 2006 Jul 6 


HD 141569 A 


2002 


Apr 21, 2008 Jul 9 


HD 150193 


2002 


Jul 22, 2007 Mar 6 


HD 163296 


2001 


Aug 6, 2001 Aug 8, 2002 Apr 21, 2002 Jul 22 


HD 190073 


2002 


Jul 22. 2005 Sep 22 


HL Tau 


2001 


Oct 25, 2002 Dec 17-18 


LkHa 330 


2002 


Dec 16, 2002 Dec 18, 2003 Nov 2-3, 2004 Dec 27, 2004 Dec 29 


MWC 480 


2001 


Jan 30, 2001 Aug 8, 2002 Jan 3, 2002 Dec 17-18 


MWC 758 


2002 


Apr 21, 2002 Dec 16, 2002 Dec 18 


SR9 


2004 


Jul 23-24 


SR 21 


2002 


Apr 21, 2006 Jul 6-7, 2007 Mar 06 


SU Aur 


2002 


Dec 16, 2002 Dec 18, 2003 Nov 2-3 


T Tau 


2002 


Jan 3, 2002 Dec 17, 2004 Dec 27, 2008 Dec 10 


TW Hya 


2002 


Dec 18-20, 2004 Dec 27, 2005 Dec, 2005 Apr 24, 2006 Apr 7-8 » 


UX Tau A 


2006 


Dec 28-29, 2007 Oct 29-30 


VV Scr 


2001 


Aug 6, 2003 Jul 9-11, 2004 Jul 23-24 


V1121 Oph 


2002 


Jul 21, 2008 Apr 18 


Wa Oph 6 


2005 


Apr 26 



* 2006 April observations obtained with Phoenix on Gemini South. 
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TABLE 2 

Stellar, Parameters 



Star 




Li, 


d 


i 


'''dust 


log(Af) V sin(j) Type^ 


Refs 




(Mq) 


(Lq) 


(pc) 


n 


(AU) 


(Mq yr~i) (km s'^) 





AA Tau 


0.67 


0.98 


140 


AB Aur 


2.4 


47 


144 


AS 205 A 


1.0 


4.0 


125 


BP Tau 


0.77 


0.83 


140 


DF Tau 


0.53 


2.97 


140 


DL Tau 


0.56 


0.77 


140 


DoAr 44 


1.4 


1.3 


125 


DO Tau 


0.72 


1.38 


140 


DR Tau 


0.4 


3.00 


140 


iau A 


0.73 


1.50 


140 


Lriv iau 


U. / 


1 1 '7 
1.1/ 


1 /I n 
14U 


OAT A « 

CjM Aur 


0.5 


0.74 


140 




0.6 


1.0 


125 




1 H 


D.o 


o4 


UPl 1/11 c^fiO A 


z.uu 


zo. / t 


1 HQ 
lUo 


HD 150193 


2.3 


1.47 


150 


HD 163296 


2.3 


36.0 


122 


HD 190073 


5.05 


470.8 


767 


HL Tau 


0.55 


0.9 


140 


LkCa 15 


0.7 


0.74 


140 


LkHa 330 


2.5 


16 


250 


MWC 480 


1.65 


11.5 


140 


MWC 758 


1.80 


11 


140 


SR 9 


1.2 


2.7 


160 


SR 21 


2.0 


11 


125 


SU Aur 


1.97 


10.70 


140 


T Tau 


2.41 


8.91 


140 


TW Hya 


0.7 


0.25 


51 


UX Tau A 


1.1 


1.0 


160 


VV Ser 


2.6 


49 


260 


V1121 Oph 


0.9 


1.5 


125 


Wa Oph 6 


0.9 


2.9 


125 



75± 10 




-8.2 


21± 


0.30±0.01 


-5.8 


25±10 


O.lSiO.Ol 


-7.1 


30I2 


0.08±0.05 


-7.9 






-6.7 


25± 5 




-7.6 


45±10 






42± 4 




-7.5 


37± 3 


0.07±0.03 


-5.1 


■^7-1- 1 




— 1 .0 


52ibl0 




—9.3 


51it 2 


n nn+0.08 


—8.2 


46it 7 




— 7.2 


01 -L 1 n 
zlitlU 




■ ■ ■ 


51ib 3 




— 11.0 


38± 9 




-6.2 


51± 2 


0.28±0.01 


-7.1 




0.62±0.01 




53± 1 




-8.8 


58± 4 


'-'■J^'-'-o.oi 


-8.8 


42±10 




-8.8 


26± 7 


0.28±0.01 




16± 4 


0.35±0.03 




34±10 


23+° " 


-7.5 


22±10 






52±10 


0.18±0.03 


-8.2 




0.22±0.05 


-7.5 


7± 1 


0.06±0.01 


-9.4 


29±10 




-9.0 


70± 5 


0.59±0.07 


-5.2 


38±10 




-7.0 


39±10 




-7.0 



11.4 


TT 


I6I13I25I26I42I 


80 


H 


I14I19I22I33I38I 


14.9 


TT 


I5I17I20I 


7.8 


TT 


I4I13I26I 


16.1 


TT 


I6I7I13I25I42I 


16 


TT 


I6I13I25I26I42I 




TT, Tr 


527 




11.1 


TT 


I6I13I26I41I42I 


< 10 


TT 


I3I4I13I23 


25| 


10 2 


TT 


1131 125137 


42| 


18.7 


TT 


I8I13I251I42I 




rprp rp 

1 1 , ir 


I4I13I23I 




rprp 
1 1 


I5I27I 






H Tr 


|9lT0l 




258 


H, Tr 


I28I31I40I 


100 


H 


I21I22I29I39I 


120 


H 


I20I23I33I 


15 


H 


I15I20I30I33I 




TT 


I6I26I41I42I 


12.5 


H, Tr 


I4I13I23I 




H, Tr 


I9I10I36I 


85 


H 


\vmo\mi'7\ 




H 


1121201 




15.2 


TT 


|7|17| 






TT, Tr 


|5|27| 




65.0 


TT 


I3I11I13I 


20.1 


TT 


I2I13I25I 


4 


TT, Tr 


I16I18I23I35I 


25.4 


TT, Tr 


I9I13I24I32I42I 


200 


H 


I19I22I34I 




TT 


527 




22.9 


TT 


I5I17I27I 



References. — (2) Akeson et aL 2002 (3) Akeson et al. 2005b (4) Akeson et al. 2005a (5) Andrews et al. 2009 (6) 
Beckwith et al. 1990 (7) Bouvier 1990 (8) Bouvier et al. 1995 (9) Brown et al. 2007 (10) Brown et al. 2009 (11) Calvet 
et al. 2004 (12) Chapillon et al. 2008 (13) Clarke & Bouvier 2000 (14) Corder et al. 2005 (15) Cuttela & Ringuelet 
1990 (16) do la Reza & Pinzon 2004 (17) Eisner et al. 2005 (18) Eisner et al. 2006 (19) Eisner et al. 2007 (20) Eisner et 
al. 2009 (21) Fukagawa et al. 2003 (22) Hillenbrand et al. 1992 (23) IsoUa et al. 2009 (24) Johns-Krull et al. 1998 (25) 
Johns-Krull & Gafford 2002 (26) Kitamura et al. 2002 (27) Lahuis et al. 2007 (28) Merm et al. 2004 (29) Monnior et al. 
2005 (30) Montcsinos et al. 2009 (31) Mora et al. 2001 (32) Najita ot al. 2007 (33) Pogodin et al. 2005 (34) Pontoppidan 
et al. 2007 (35) Qi et al. 2004 (36) Salyk ot al. 2009 (37) Simon ot al. 2000 (38) Tannirkulam et al. 2008 (39) van don 
Ancker et al. 1998 (40) Weinberger et al. 1999 (41) White & Ghez 2001 (42) Wichmann et al. 1998 
^ Classifieation used in this work. H:Herbig Ac/Be, TT:T Tauri and Tr:Transitional. 



have strong central depressions that are hkely inconsis- 
tent with simple disk emission from a Keplerian disk, 
with the latter, unfortunately, having no observable lines 
with Jup >20. Our sample also includes a few sources 
known to be strongly centrally peaked (AS 205 N and 
DR Tau; Bast et al. 2011), with the overall line shape 
consistent with the sum of a Keplerian profile and a 
slow disk wind (jPontoppidan et all 1201 ID . We do not 
make an effort to model or understand these exact struc- 
tures, and instead focus here on a derivation of inner 
radii from the line wings. For more in-depth studies of 
overall line shapes, utilizing higher-r esolution data from 
VLT-CRIRES, we d irect t he reader to lPontoppidan et alJ 
(|2ml) . lBairetn] (poTTI ). and J.M. Brown et al. (2011, 
in preparation). 

4.2. Profile Modeling 

4.2.1. Procedure 

While temperature profiles have been measu red for 
outer disk dust (e.g., Andrews et al- ,2009; .Isella et al.l 



120091) ■ the temperature profile for inner disks or for 
molecular line-emitting layers is unmeasured. In fact, 
it is believed that the disk upper atmosphere is ther- 
mally decoupled from the dust, and its temperature is 
set by a complex balance of gas heating and coolin g 
(IGlassgol d fc Naiital [200l IKamp fc Dullemondl [200l . 
Therefore, we have chosen a model that makes a min- 
imal number of assumptions about the underlying tem- 
perature. In addition, our model was designed to provide 
a robust way to determine the CO inner radius, without 
the measurement of this radius depending strongly on 
the choice of model. In other words, it allows for a good 
fit to the line wings (and hence i?in) independently of 
lower- velocity portions of the line profile. 

We have chosen to model the line profiles as emission 
from a disk with CO luminosity per unit radius Lco{R), 

where Lco.entire disk = Lco{R)dR and Lco{R) is 

a broken power law with BP for i^in < R < Rmid and R'^ 
for i?mid < R < Rout- By directly fitting for the shape 
of Lco(-R)) rather than beginning with the temperature 
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TABLE 3 
Line Profile Parameters 



Star 


Lines in Composite 


FWHM (km s^^) 


-^^gauss (AU)^ 


AA Tau 


P 


[30,31,32,37,38,39,40) 


92 


0.09 


AB Aur 


P 


[26,27,28,29,30,31,32,34,36,37,38,39) 


28 


0.59 


AS 2(15 N 


P 


[22,30,31,32,33,34,36,37,38,39,40,41,42,43,45) 


43 


0.12 


BP Tau 


P 


[36,37,38) 


87 


0.03 


DF Tau 


p 


[30,31,32,36,37,38,40) 


79 


0.10 


DT, Tail 


p 


[30,31,32) 


87 


0.02 


DoAr 44 


p 


30,31,32,34,35,36,37,38,40) 


61 


0.24 


DO Tail 


p 


[30,31,32,33,36,38,39,40) 


87 


0.05 


DR Tau 


P 


[30,31,32,33,34,36,37,38,39,40) 


29 


0.27 


GG Tau A 


p 


30,31,32,33,34,36,37,38,39,40) 


68 


0.07 


GK Tan 


p 


[30,31,32,33,36,37,38,40) 


97 


0.06 


GM Aur 


p 


[9,10,11,12,14) 


47 


0.19 


GSS '^Q 


P 


[30,31,32,37,38,40) 


80 


0.06 




p 


[1,2,3,5,6,7,8,9,10,11,12) 


19 


1.49 


HD 1 41 5fiq A 


p 


[1,2,3,4,5,6,7,8),R(1,0) 


24 


3.68 


HD 150193 


P 


^30, 32, 36, 37) 


53 


0.40 


Hn 1 fi32Qfi 


P 


[27,30,31,32,37,38) 


83 


0.26 


HD 1 9007'^ 


P 


[30,31,32,34,36,37,38) 


24 


3.38 


HL Tau 


P 


[27,30,31,32,33,34,36,37,38,39,40) 


96 


0.05 


LkHa 330 


P 


[30,31,32,33) 


28 


2.10 


MWC 480 


P 


[26,27,28,29,30,31,32,36,37,38,40) 


72 


0.08 


MWC 758 


P 


[30,31,32,33,36,37,38) 


32 


0.19 


SR 9 


P 


[30,31,32,37) 


77 


0.08 


SR 21 


P 


[5,6,7,8,9,10,11,12) 


18 


1.98 


SU Aur 


P 


[30,31,32) 


121 


0.10 


T Tau 


P 


[30,31,32,33,34,36,37,38,39) 


62 


0.19 


TW Hya 


P 


[6,7,8,9) 


17 


0.05 


UX Tau A 


P 


[1,2,3,4,5,6,7,8,10,11),R(1,0) 


21 


1.04 


VV Ser 


P 


[30,31,32,34,36,37) 


64 


0.72 


V1121 Oph 


P 


[31,32,36) 


109 


0.04 


Wa Oph 6 


P 


[30,38) 


139 


0.02 



^ Derived from the velocity at 1.7 X HWHM and disk inelination. See the text for details. 

star and the upper limit producing maximum velocities 
similar to the NIRSPEC resolution. We also tested a 
similar range of values for ii!mid, though i?mid was by 
definition always larger than or equal to i?in. 

4.2.2. Model Results 

Most of the model parameters turn out to be strongly 
degenerate, so that they cannot be uniquely determined 
from the line profiles, and since our focus was to study 
CO inner radii, we did not make an effort to understand 
these parameters in detail. Instead, our goal was to in- 
vestigate whether i?in could be determined robustly. We 
find that if p is strongly negative, then the line shape be- 
comes very sensitive to Rm- This is demonstrated in the 
middle panels of Figure [21 which show as a function 
of and p, with all other parameters being free vari- 
ables. If p < —1.5, then i?in is robustly determined. If p 
is (i.e., LcoiR) is fiat from Rin to i?mid), then the line 
shape becomes insensitive to Rin, but is instead sensitive 
to i?mid- (For intermediate values of p, the line shape is 
sensitive to a combination of i?in and i?inid-) Thus, the 
line profile is sensitive to steep drops in Lco(-R)- This 
is due to the fact that in order to be sensitive to some 
i?o (either R^ or i?mid), any change to this parameter, 
ARq, must result in a significant change to the total line 
luminosity, i.e., L(Rq)ARo must be a significant fraction 

of JI^""^* LcoiR) dR. Figure [21 we show line profiles 
with p = and p = —1.5 for two different values of i?in, 
and demonstrate that with p = —1.5 the line profile is 
very sensitive to i?in, but with p = 0, the line profile 



and CO column density profiles, T(i?) and N{R), and 
using these to deduce Lqo (^)j we make the problem both 
computationally manageable and simple to int erpret. A 
similar procedure was used by (jCarr et al.ll2Q04f ) to model 
CO and H9O emission from the disk around the young 
star SVS 13. 

A broken power law is empirically convenient, as it has 
enough free parameters to properly fit most line profiles 
(in many cases a single power law is not sufficient to 
fit the observed line profiles), and yet the inner radius 
turns out to be well constrained. A broken power-law 
approximation for Lco(-R) also has a physical basis. At 
aU radu, LcoiR) ~ ^iR) x 27ri?, where F{R) is the fiux 
emitted by the CO per unit radius. At small radii, if 
the gas is optically thick, temperatures are high and so 
the emission near 5 fiia is in the Rayleigh- Jeans regime 
(i.e., F cx T). At larger radii, however, there should be 
a steeper decrease in FiR) as temperatures drop, the 
blackbody curve peaks at or beyond 5 /im, and F drops 
much more rapidly with T. 

In all models, we have assumed i?out = 100 AU, as 
the models are insensitive to i?out beyond ~ a few AU. 
We allow q to take values of —1.5, —2, —2.5, or —3. 
q primarily affects the degree of double peak in the line 
profile, but is not well constrained because of degeneracy 
with other parameters. We tested values for p from —2 to 
and, as we will discuss, this parameter has an important 
effect on the determination of the CO inner radius, Ri^. 
Values of i?in from 0.005 to 10 AU were tested, with the 
lower limit being the approximate radius of a solar mass 
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Fl(i. 2. Left: composite emission lines and best-fit models witfi p = —1.5. Middle; contours for best models as a function of Rin 
and p. witli dark representing lower ■ (Other model parameters are left as free variables). Squares show the best-fit R^^ for p = —1.5; 
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TABLE 4 
Model Fits 
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AA Tau 


0.10 


0.20 


—3.0 


0.010 


0.20 


—3.0 


AJd Aur 


0.40 


0.50 


— 1.5 


0.050 


1.00 


— 1.5 


A c one "XT 


0.07 


0.10 
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0.010 


0.10 
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BP Tau 


0.03 


2.00 
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0.20 
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DF Tau 
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Fig. 3. — Normalized line profiles (solid lines) with p = (left) 
and p = —1.5 (right) for two different values of 0.005 AU 

(black) and 0.1 AU (gray). All models assume q = —2, -Rmid = 0.5, 
R^^t = 100 AU, Mi, = Mq and i = 45°. Dashed and dotted 
lines show the contributions from to Rmid a^nd -Rmid to Rout, 
respectively. Note that with p = the line profile is relatively 
insensitive to while with p = —1.5 it is very sensitive to iJin. 

changes only slightly with changes in Rin- 

Due to the dependence of and i?niid on p, we have 
explored two classes of solutions that can fit the data 
well in most cases. In the first approach, p is set to —1.5, 
and Rin is well-determined. In the second p is set to 0, 
and -Rmid is well-determined. The best-fit parameters in 
these two cases are presented in Table H) The best-fit line 
profiles and corresponding Lcq{R) profiles are shown in 
Figure [H The p = case places no constraints on the 
inner edge of the CO emission and has a sudden steep 
drop-off in Leo (R) occurring at i?mid ■ The p = —1.5 case 
has a distinct region free of CO between i?,^ and Rin- We 



o 

II 




.001 



100 



Fig. 4. — Plot of Rin with p fixed at —1.5 and -Rmid with p fixed 
at 0. -Rmid with p = is typically 1 — 5 X .Rin with p = —1.5 (see 
dashed lines). 

believe that the second class of models is more physically 
motivated, as CO can be cleared out near the star due 
to photodissociation or disk truncation. However, we 
cannot empirically rule out the other class of models. 
There are a few sources in which the p = solution 
results in a significantly-improved fit — most notably AS 
205 N and Wa Oph 6 (see Figure E]). However, AS 205 
N is known to be a membe r of a small subs et of disks 
with "peaky" line profiles ()Bast et al.l 120 lit) , in which 
low- v elocity flux may be enhan ced by a low-velocity disk 
wind (iPontoDpidan et a"Lll2011|) . Thus, a Keplerian disk 
model is probably not appropriate in this case, and we 
do not believe that the model with p = is necessarily 
more physically realistic than the model with p = —1.5. 

In most cases, one class of fits is not preferred over the 
other, and since we feel the p = —1.5 case is more phys- 
ically motivated, we adopt this model for the remainder 
of this work. However, in Figure |4l we show that -Rin 
with p — —1.5 and -Rmid with p = are related such 
that -Rmid ^ 1 — 5 X -Rin (because with p — there is 
still some fiux contribution from the region inside -Rmid)- 
Therefore, uncertainty about the choice of model simply 
means that there is a systematic uncertainty in the loca- 
tion of the drop-off in Lqq by a factor of a few, and any 
trends in -Rin should be robust. 

When comparing -Rin with other disk radii elsewhere 
in the paper, we will refer to our best-fit -Rin as -Rco- 

4.2.3. Error Estimates 

Assuming zero error in disk inclination, it is apparent 
from the diagrams in Figure [2] that i?in is uncertain 
to ^ 0.5 dex, or factors of a few. Uncertainties are also 
introduced due to errors in disk inclination. These scale 
roughly as AR/R ~ cosiAi, and so are largest for low- 
inclination disks. With a typical inclination (45°) and 
error (±10°), errors in -R (derived from Kepler's law) are 
-1-50% and —25%. Thus, they are similar to or somewhat 
smaller than the systematic uncertainties. However, for 
smaller inclinations, uncertainties in i may dominate. In 
any case, an uncertainty of 0.5 dex should reasonably 
be assumed for .Rin. 



10 



Salyk et al. 



For a few sources with narrow line profiles, it is appar- 
ent from the lack of a large radius cutoff in ^ (Figure 
[2]) that the models are fairly insensitive to differences in 
i?in at radii larger than the best-fit value. These include 
AB Aur, HD 135344 B, HD 141569 A, HD 190073, LkHa 
330 and SR 21. 

As discussed in Section l4.2.21 there is a systematic un- 
certainty introduced by the two classes of models capable 
of reproducing the observed line shapes. There is also a 
potential bias introduced by the choice of lines used in the 
composite line profile. In theory, one should be able to 
obtain the same i?in from all emission lines and we have 
tested whether using only low-excitation lines (J < 12) 
yields the same result. In practice, we find that utilizing 
low-excitation lines only, the derived i?in are similar to 
or larger than the nominal i?in, with the ratio typically 
being 1-2. Although biased towards larger i?in, this is 
within our known uncertainty of a factor of a few, and 
we believe that the high-excitation lines are likely to be 
more sensitive to i?in, and so yield the more accurate fit. 

4.2.4. Benchmarking to Other Work 

At least two other methods for determining i?in have 
appeared in the literature. In one, a CO temperature 
(and density) profile is assumed, and Rm is fit as a sin- 
gle free variable in a d i sk model (e.g.. IBlake fc BoogertI 
I200I ISalvk erall[2009l: IBast et al.l[201lir In our model, 
in contrast, the temperature, density, level populations, 
and emitting area are all wrapped up into Lco{R)- To 
compare the two models, we fit model line profiles, con- 
struct ed using the procedure described in ISalvk et"al] 
(|2009| ). with the simple model described here, setting 
p = q. We find that disk temperature profiles of the 
form T oc Rp with p = —0.6 are equivalent to models 
with Lco{R) oc R^^, and T oc is equivalent to 
~Lco(i?)oci?-i-5. 

Compared to the values derived bv ISalvk et al.l ()2009l ). 
we find similar radii for many sources, with the notable 
exception of HD 1 41569 A, LkHa 330, and SR 21, for 
which Sa lvk et al.l (2009) find radii larger by factors of 
a few to ^ 10. All three of these sources, however, are 
ones in which is not sensitive to Rm at radii larger 
than the best-fit valu e. We also find radii consistent 
with IBast et"all (|20T1 for AS 205 A and TW Hya, who 
estimate inner radii of 0.04 and 0.1 AU, respectively (al- 
though they caution that the parameter space for these 
models was not well explored; J. E. Bast, private com- 
munication). Our respective results are inconsistent for 
VV Ser, for which Bast et al. derive an inner radius of 
0.08 AU and we derive an inner radius of 0.72 AU; how- 
ever, a closer look at their model reveals that they are in 
a regime similar to our p = case, in which is simply 
not very sensitive to the choice of inner radius. 

Another common approach to estimating Rin is to sim- 
ply choose some velocity (typically either 2x the half- 
width at half maximum (HWHM), or the half- width at 
zero intensity), and set Rin to the radius with that Ke- 
plerian velocity. The difficulty of this approach is that 
it is not obvious which velocity to choose. In Figure [5l 
we compare i?in to inner radii derived from fitting the 
profiles with single or double (emission plus absorption) 
Gaussians (i?gauss)- In particular, we show solutions in 
which i?in is cierived from the velocity at 1.7x HWHM, 
which comes closest to reproducing our results. Devia- 



tions from 1:1 are of order a factor of a few, and so this 
simple approach is remarkably consistent with our more 
complex model. Thus, we suggest that using the velocity 
at 1.7x HWHM is a reasonable choice for calculating _Ri„ 
using simple Gaussian fits. We use this result to incor- 
porate LkCa 15 into our a nalys is, utilizing the FWHM 
measured by iNaiita et al.l ()2003() . This result also gives 
us confidence that the model-derived Rm is reflecting the 
bulk line shape and FWHM, and is not instead some spu- 
rious result heavily biased by the line/ continuum ratio, 
the noise level, or any other aspects of the data. 

4.2.5. Comparison with Spectro-astrometric Results 

Four sources in our sample (HD 135344 B, SR 21, 
TW Hya, and VV Ser) have inner radii derived from a 
combined line shape and spectro-a s trome tric (SA) profile 
analysis (Pont oppidan et al.|[2008l 1201 ID . Ahhough the 
overlapping sample is small, and two of these are sources 
for which we do not have a strong upper limit on R^ 
(see Section l4.2.4|) . our results ap pear broadly consistent 
with these results. Pontoppidan et aD (|2008[ ) measured 
CO inner radii for HD 135344 B , SR 21 and TW Hya 
assuming a power-law disk temperature profile, and de- 
rived values within a factor of 3.5 of our Rin- Note that 
the SA signal is most sensitive to the extent of the CO 
emission, and also provides a hard upper limit on 
However, the SA profile is not very sensitive to Rm, and 
so modeling that incorporates SA is subject to the same 
uncertainties in the gas temperature profile, and does not 
n ecessarily determine Rm m ore accurately. 

iPontoppidan et all ()2011[ ) report SA radii, which they 
define as the radius at the peak of the SA profile, for 
HD 135344 B, SR 21, TW Hya and VV Ser. The SA 
signal is the flux-weighted mean position of the emission 
at each velocity and is therefore sensitive to the distribu- 
tion of emission rather than simply the inner boundary. 
The SA radii will therefore always be larger than R^ 
unless the emitting region is infinitesimally thin. In the 
context of our two-power-law model, the SA radii would 
be affected by p, q and i^mid, since these all affect the 
amount of flux at large radii. Thus, the SA profiles are 
complementary to the results from line shape analysis, 
and the two can potentially be used in concert to derive 
the shape of Lc,o(R). T he SA radii (Rsa) derived by 
([Pontoppidan et aD 120111 ) are factors of 2-5 larger than 
our i?in . We have tested our two-power-law models with 
a simple code to calculate Rsa and find that we can si- 
multaneously reproduce the observed Rm and Rsa by ad- 
justing other model parameters. Therefore, our Rm are 
consistent with and complementary to the SA results. 

5. LINE FLUX MODELS 

Rotation diagrams can be used to derive characteristic 
column densities, emitting areas and temperatures for 
the CO emitting layer. Rotational levels are assumed 
to be populated according to LTE. Although these pa- 
rameters are not fully realistic, since the emission actu- 
ally comes from a range of radii and heights in the disk 
atmosphere, they provide a convenient way to roughly 
characterize and compare the emission within a sample 
of disks. In this work, we will also compare these param- 
eters with i?co- 

5.1. Fitting Procedure 
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Fig. 5. — Modeled inner radii compared with inner radii derived 
from Gaussian fits and velocities at l.TxHWHM (with the factor 
of 1.7 determined empirically). 

We fit CO line fluxes with an LTE slab model fo l lowin g 
the procedure described in detail in iSalvk et alj (|2009f ). 
Sources with fewer than two measurable line fluxes with 
Jup > 20 or a significant amount of self-absorption in 
low-excitation lines were excluded from the analysis, as 
these did not provide reliable model fits. Line fluxes 
were calculated using Gaussian fits to the emission lines, 
and rotation diagrams were fit using a grid of LTE slab 
models with a single temperature (T), CO column den- 
sity {N), and emitting area (A) as free parameters. In 
addition, ^^CO detections or non-detections were used 
to constrain N and remove model degeneracies. Best-fit 
model parameters are listed in Table [5] 

6. DISCUSSION 

6.1. Correlation with Corotation Radii 

For a sample o f five T Tauri stars observed by 
INaiita et all jSool), CO inner radii were found to be in 
the range ~ 0.5 — 1 times the corotation radius, the ra- 
dius at which the disk's Keplerian angular velocity equals 
that of the stellar surface: 

_ GAhRlsin^ii) 
^ ~ (i;sin(z))2 ■ 

For convenience, we assume that i^, = idisk, although 
the two need not be exactly the same. Rc is expected 
to be similar to the magnetic truncation radius, within 
which the disk is cleared by magnetospheric accretion 
()Shu et al.lll994l ). In Figure El we plot Rco against i?c, 
as well as lines representing Rco = Rc and Rco = 
0.5 X Rc- Symbol sizes are proportional to the stellar 
mass. We find that Rco is consistent with being ~ 0.5- 
1 X the corotation radius for T Tauri disks, but that this 
relationship quickly breaks down for transitional disks 
and disks around HAeBe stars. Additionally, there is 
a hint that the degree of discrepancy scales with stellar 
mass, a trend we discuss in Section [6.21 



TABLE 5 
Rotation Diagram Fit Results 



Star 






logM (AU2)) 


AA Tau 


950 


18.6 


-1.3 


AB Aur 


600 


18.8 


0.3 


AS 205 N 


975 


18.7 


-0.3 


DF Tau 


1025 


19.1 


-0.7 


DL Tau 


1675 


18.6 


-1.7 


DoAr 44 


1150 


18.2 


-1.6 


DO Tau 


1575 


19.0 


-1.7 


DR Tau 


1250 


18.9 


-1.0 


GG Tau 


1500 


18.6 


-2.1 


GK Tau 


1600 


18.1 


-1.5 


GSS 39 


1675 


18.1 


-1.4 


HD 135344 B 


900 


17.9 


-1.5 


HD 141569 A 


275 


18.3 


2.0 


HD 150193 


700 


18.1 


-0.5 


HD 163296 


825 


18.5 


-0.5 


LkHa 330 


850 


17.4 


-0.6 


MWC 480 


975 


18.0 


-0.8 


MWC 758 


700 


18.3 


-0.2 


SR 9 


1575 


16.7 


-0.8 


T Tau 


1325 


17.9 


-0.6 


TW Hya 


700 


17.7 


-1.9 


V1121 Oph 


1125 


18.5 


-1.4 


Wa Oph 6 


1675 


18.3 


-1.6 



^ It is not straightforward to define error bars, because of 
the degeneracy of model parameters iSalvk et al.lf2b09t) . Typ- 
ical uncertainties are ±200 K for T, and ±0.5 for log(A'^) and 
log(A). 
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Fig. 6. — CO inner radius against corotation radius. The size of 
the symbol is proportional to the stellar mass. The solid line shows 
a 1:1 correspondence, while the dotted line has Rco = 0.5xiJcorot. 



There exists a known correlation between inner dust 
disk size and stellar luminosity, in which inner dust disk 
sizes appear to be determi ned by dust sublimation, with 
Tsnh ~ 1000 - 1500K fsee iDullemond fc MonnieillMol . 
and references therein). For blackbody grains in a disk 
with an optically thin inner region. 



6.2. A Size- Luminosity Relationship for Gas Disks 



(2) 
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This formula neglects many complications, including 
grain properties, the pressure dependence of sublimation 
temperatu re, and grain scattering , whic h are discussed 
in detail in lDuUemond fc Monnieif ()2010t ): since the rela- 
tive importance of these various effects has not yet been 
resolved, we work with this simple version. 

CO inner radii might also be set by dust sublimation, 
since CO can be dissociated by UV radiation if grains 
are not available for shielding. However, previous stud- 
ies of CO inne r radii were lim ited to small ranges in stel- 
lar luminosi ty (iNaiita et al.f 2003: Blake & Boogert 2004; 
ISalvk et al.ll2009D . In addition, low-mass cTTs's have 
dust sublimation radii that are similar in size to the stel- 
lar corotation radii, so distinguishing between the two 
effects is difficult. In Figure [71 we show the relationship 
between luminosity and -Rco for our complete sample of 
disks. Note that instead of simply plotting the stellar 
luminosity, we also include the accretion luminosity if 
the accretion rate is known, as this is the total luminos- 
ity seen by the disk. W e assume Lace = O.SGM^^Af/i?* 
(jHartmann et al.lll998t ) . Neglecting to include the accre- 
tion contribution to the luminosity yields outliers with 
high accretion rates (see Figure [5]). 

In Figure [71 we show the strong correlation between 
Rco and L^, + Lace, consistent with dust sublimation 
temperatures of T - 1500 - 2000 for cTTs disks and 
slightly lower temperatures, ^ 1000 — 1500 for HAcBe 
disks. Transitional disk radii are larger, but also follow a 
similar trend (discussed further in Section |6.4[) . Includ- 
ing all disks, the p- value associated with linear regression 
of log(i?co) against log(L^, -I- Lace) is 6 x 10^^, and ex- 
cluding transitional disks, it is 10~^°, so the increase in 
radius with luminosity is highly statistically significant. 
Therefore, we believe this is strong evidence for a depen- 
dence of Rco on dust sublimation. 

The slope of the trend (excluding transitional disks), 
0.7 ± 0.1, is somewhat steeper than the slope of 0.5 ex- 
pected from dust sublimation alone, at the 2 a level. 
In more detail, Rco may be set by a balance of pho- 
todissociation (which eats outward in the disk) and ac- 
cretion (which replenishes inward). This could result in 
relatively smaller Rco for T Tauri disks and larger for 
HAeBe disks due to the respectively lower and higher 
photodissociating UV fiuxes. 

6.3. Comparison with Interferometric Dust Inner Radii 

If dust acts as a shield against photodissociation, then 
we should also expect a correlation between Rco and 
measured dust inner radii from near-IR interferometers. 
We show these radii in Figure [9l along with a line mark- 
ing a 1:1 correlation, and another marking the best linear 
fit. Although there is a statistically significant correla- 
tion between the two variables, there is also significant 
scatter, of order 0.5 dex. Uncertainties in i?co can be 
a factor of a few (see Section I4.2.3P , which may account 
for some of the scatter. Also, dust inner radii are mea- 
sured using different models in different studies (includ- 
ing thin rings, disks, both inclined and not inclined), and 
so it is possible that some of the scatter arises from the 
choice of dust model, or from dust models not accounting 
for the disk inclination. Using a coherent sample of in- 
terferometric visibilities, analyzed in the same way, and 
accounting for disk inclination, could test this hypothe- 



sis. Another possibility is that there is real scatter due 
to different rates of photodissociation and/or replenish- 
ment via accretion. This is consistent with the obser- 
vation that several HAeBe disks have anomalously large 
i?cO: while several cTTs's have anomalously small -Rco- 

6.4. Transitional Disks 

As discussed in several prior st udies (I Naiita et al.l 
l2003HRettig et al.|[200llSalvk et al.| [2009). Rco lies weh 
inside -Rtrans for many transitional disks (where -Rtrans 
is the radius at which the disk becomes optically thick). 
This result has also been confirmed yia spectroastrome- 
try fo r HD 135344 B and TW Hya (jPontoppidan et al.l 
20081). In F i gure [lOl we closely reproduce the results of 
Salvk et"al] ( 20091) . who showed that the best- fit -Rco lie 
somewhere between -Rsub and -Rtrans for all disks. An 
important caveat is that, as discussed in Section 14.2.21 
-RcOi the upper limit to -Rco is not well constrained for 
HD 135344 B, HD 141569 A and SR 21. The emission 
in SR 21 and HD 141569 A could originate at or near 
- Rt.ra.ns, as is obscrvcd for the transitional disk HD 100546 
(|Brittain et al.ll2009f ): however, for HD 135344 B, some 
of the CO gas mus t originate from within th e inner, op- 
tically thin region (jPontoppidan et al.l[2008l). 

Using the procedures described in Section li. 2. 1[ we de- 
rive significantly larger radii for transitional disks than 
for classical disks at the same luminosity. This difference 
is typically an order of magnitude — much larger than 
our uncertainties — and holds true for almost three or- 
ders of magnitude in luminosity. Here we discuss some 
possible explanations for this difference. 

One possible explanation is that this is the result of 
some systematic bias in our analysis. Two possibilities 
come to mind. First, since -Rco depends on the line 
wings, one might derive smaller inner radii by "fitting 
the noise" at large velocities. However, since transitional 
disks tend to have lower continuum S/N level s, and lower 
line / continuum ratios than classical disks (jSalvk et al.l 
120 111) , this would tend to bias one toward deriving larger 
characteristic velocities and hence smaller inner radii for 
transitional disks. In addition, we find that the transi- 
tional disk radius discrepancy holds even if we examine 
-Rin derived from the HWHM. A second possibility is that 
-Rco appears to depend somewhat on the choice of lines 
included in the line composite, with lower-excitation lines 
yielding larger radii by up to a factor of ^ 2. Since 6/8 
transitional disks in our survey (GM Aur, HD 135344 B, 
HD 141569 A, SR 21, TW Hya, and UX Tau A) did not 
have high-excitation emission lines, their low-excitation 
lines were analyzed instead, and this could result in spu- 
riously large radii for these disks. However, this bias 
produces only up to a factor of two difference in -Rco, 
while the observed difference is an order of magnitude. 
Furthermore, the discrepancy is also seen for DoAr 44 
and LkHa 330, which do have high-excitation emission 
lines. Instead, we suggest that the lack of high-excitation 
emission is a reflection of the relatively larger emitting 
radii (and thus lower emitting temperatures) for these 
disks. 

A physically-motivated explanation for the radius dis- 
crepancy is that -Rco in transitional disks is not set by 
sublimation, but rather by dynamical truncation by an 
embedded protoplanet. Although embedded protoplan- 
ets have been posited as a possible explanation for the 
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Fig. 7. — CO inner radius against luminosity. Diamonds are cTTs disks, squares are HAeBe disks and stars are transitional disks. 
Filled symbols have measured accretion rates, which arc incorporated into Lacci while unfilled symbols have Lace = 0. Dashed lines show 
theoretical curves for dust sublimation radius as a function of luminosity. 
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inner clearings in transitional disks, they are expected 
to orbit at radii capable of producing the sharp transi- 
tion between the optically thin inner and optically thick 
outer disk regions. Since the Rco we find here is usually 
significantly smaller than i?trans we would actually need 
to invoke dynamical truncation by additional planets, in 
smaller orbits. In particular, we w ould predict planets lo- 
cated at radii near ~ 0.5 x Rco (|Artvmowicz fc Lubowl 
I1994D . This is in line with recent work suggesting that 
multiplanet systems may be required to explain transi- 
tional disks (|Zhu et al..,2011; .Dodson-Robinson fc Salvk 



1201 If ). However, this explanation would need to be con- 
sistent with the apparent increase in i?co with luminos- 
ity. This is not unreasonable, as the mass and location of 
the protoplanets may depend on the disk and/or stellar 
mass. 

Another possible explanation for the radius discrep- 
ancy, and one qualitatively consistent with the increase 
of Rco with luminosity, could be that the dust disk is 
still truncated at i?sub, but that photodissociating UV 
photons penetrate a finite distance into the tenuous tran- 
sitional dust disk, making Rco > Rsuh- A possible quan- 
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titative test for this hypothesis would be to calculate a 
UV penetration depth and compare to the radius dis- 
crepancy. However, this is difficult to test in practice, 
since the local dust density at i?sub is not known. (Al- 
though SED models can place some constraints on the 
location and surface density of the inner disk dust, they 
do not do so to the precision required here.) Yet, with 
sufhciently low dust densities, this explanation is plau- 
sible; for example, with k ^ 10^ cm^ and a dust 
density p = 10~^® g cm~'^ (equivalent to a vertical dust 
surface density of ~ 10~^ g cm"'^ for a scale height of 0.1 
AU), the pe netration depth (to r = 1) is ~ 7 AU . For 
comparison, [Eisner. Chiang, fc Hillenbrandl (j2006[ ) mea- 
sure a dust surface density of 6.3 x 10~^ g cm~^ at the 
inner edge of the TW Hya disk. 

A final possibility we consider is that the CO tem- 
perature structure in the emitting layer is different for 
classical and transitional disks. Since we cannot distin- 
guish between models with steeply declining line lumi- 
nosity, Lco(-R): beginning at Rco, and models with flat 
Lco{R) at smaller radii (see the discussion in Section 
I4.2.2[) , it is possible that some aspects of the latter model 
are more appropriate for transitional disks. This possi- 
bility could be tested with radiative transfer code that 
properly treats the gas heating and cooling as well as the 
line radiative transfer. 

It is interesting to note that 'Pontoppida n et al.[ (|2011[ ). 
in contrast, find no significant difference between SA 
radii for transitional and classical disks. As discussed 
in Section [4.2.5[ the SA radii are complementary to Rm, 
as they represent the flux- weighted inean ra dius. There- 
fore, the results of [Pontoppidan et al.l (|2011[ ) suggest that 
differences in CO emission between transitional and clas- 
sical disks are erased at large radii, and the outer extent 
of the emission is similar for the two classes of disks. 
If the outer extent is the same, but the inner radius is 
larger for transitional disks, the emission from these disks 
should come from a smaller range of radii. 
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Fig. 11. — Plot of best-fit model area against CO inner radius. 
Symbols are the same as in Figure [7| The solid line shows a 
linear fit to non-transitional sources; dashed lines correspond to 
A = 2nRAR with AR = 0.15 (top) and Ai? = 0.01 (bottom). 
(HD 141569 A, not shown, has an anomalously high, and probably 
incorrect, emitting area.) 

6.5. Inner Radii and Results from Line Flux Models 

If the CO line flux is dominated by emission from near 
Rco, we should expect a positive correlation between 
i?co and the characteristic emitting area, A, derived 
from the rotation diagram flts. In fact, as shown in Fig- 
ure [11] we flnd that for classical disks, A is approximately 
oc Rco ■ If the emission comes from a ring located at i?co 
and of thickness Ai?, then the rotation diagram results 
are consistent with Ai? — 0.15 AU. Interestingly, the re- 
lationship is quite different for transitional disks. Exclud- 
ing HD 141569 A, which has an anomalously large area 
(probably due to poorly constrained high- J line fluxes) a 
similar calculation yields Ai? = 0.01 AU. This is at least 
qualitatively consistent with the fact that the emission 
from transitional disks may come from a smaller range 
of radii, as suggested by i?in and i?sA- 

With a rotational temperature. Trot, and a location, 
Rco J we can explore the relationship between dust and 
gas temperatures. In Figure I12[ we show T^ot and 
Toff, defined as the effective temperature for a black- 
body grain fully exposed to stellar plus accretion radi- 
ation. We find no statistically significant correlation be- 
tween these two variables. In addition, while models pre- 
dict gas temperatures in the upper layers of disk atmo- 
spheres to be in excess of the thermal dust temperatures 
(jClassgold & Najita 2001; Kamp & Dullcmond 2004), 
we find that Trot is similar to, or often less than, Tcff . Al- 
though we have assumed blackbody grains, a correction 
for the relative absorption and emission efficiencies for 
disk grains would make dust tempe ratures even higher 
(e.g.. lMonnier k Millan-Gabetl[20M) . The similarity be- 
tween Trot and Tcff is also curious in light of recent results 
highlighting the need for a high gas temperature in or- 
der to produce observed H2O emission lines, which have 
critical densities simil ar to the CO lines observed here 
(|Meiierink et al.|[2OO90 . A possible explanation is that 
Trot represents a characteristic temperature for the emit- 
ting region, which is lower than the temperature at the 
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inner rim. 

We also find lower Tj-ot for transitional and HAeBe 
disks than for cTTs disks, on average. This result is 
consistent with the slightly larger values of i?co/^sub 
observed for HAeBe disks, seen in Figure [T] 

7. CONCLUSIONS 



In conclusion, we have used observations of CO emis- 
sion from a large sample of cTTs, HAeBe and transi- 
tional disks to probe inner disk structure. Using a con- 
ceptually simple parameterization of the CO emission as 
a function of disk radius, we find that CO inner radii 
can be robustly determined. These inner radii are simi- 
lar to the stellar corotation radii for low-mass stars, but 
significantly larger than stellar corotation radii for disks 
around HAeBe stars. We find a strong size-luminosity 
relationship for CO inner radii, and suggest that CO 
emission is truncated at the dust sublimation radius. 
Transitional disks are obvious outliers, with CO inner 
radii typically an order of magnitude larger than classical 
disks at the same luminosity, even though the emission 
typically arises from well within the transition radius at 
which the disk becomes optically thick. 

We also compare CO inner radii with gas character- 
istics derived from rotation diagram fits. We find that 
classical and transitional disk line fltrxes are separately 
consistent with emission from a single temperature ring 
of width 0.15 and 0.01 AU, respectively, located at the 
CO inner radius. We also find systematically lower ro- 
tational temperatures for transitional disks and disks 
around HAeBe stars, than for those around T Tauri 
stars, which is consistent with observed differences in 
the ratio of CO inner radius to dust sublimation radius. 
Finally, we find rotational temperatures similar to, or 
slightly lower than, the expected temperature of black- 
body grains located at the CO inner radius. 

The data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the 
University of California and the National Aeronautics 
and Space Administration. The Observatory was made 
possible by the generous financial support of the W.M. 
Keck Foundation. 
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